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I. IMTRODUCTIOM 

Sinca 1979 thara has been a number of papars on avaluation of energy 
release rata for thermoelasticity and corresponding J-intagral. Two sain 
approaches are developed to treat energy release rate in elasticity. The 
first is based on direct calculation of potential energy rate with respect to 
(rack length (11. The second Makes use of Lagrangian forma lisa (Mather's 
theorem) (2). The direct method was extended to the problems of 
thermoelasticity by Wilson and Yu (3), K. Kishimoto, S. Aoki, and M. Sakata 
(41 and Curtin (5]. A useful discussion on the extension of the approach to 
theraoelasticity was presented by McCartney (61. 

The second approach was elaborated by Atkinson and Saelser (7] and 
discussed in (81. In elasticity the energy release rate is expressed as a 
path independent Integral. The path lnvariancy of the similar integrals for 
thermoelastlcity is discussed in aforementioned papers. 

In [3, 4] the J-integral consist of two parts: first is a common for 

elastic problems integral along a path and the second Includes a nonvanishing 
integral over the volume (area for two dimensional applications). 

Path independent integral introduced in (5) has limitations in path 
lnvariancy as was noticed in (7): the integration paths should start and end 
at certain points. It is not discussed in (5] whether the Integral represents 
any energy release rate. 
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Another path independent integral expreeeed in teres of energy-momentum 
tensor has been derived in (7). The derivation follows the Lagrangian 
formalise of [2] and essentially depends on Lagrangian chosen. A comparison 
of this integral and that obtained in present paper (see also (6J) indicates 
that an entropy tore is omitted. 

In this paper we consider translational and expansional energy release 
rates employing the formalism of irreversible thermodynamics and Crack Layer 
Approach. A damage zone usually precedes and surrounds a propagating crack. 

A system of a crack and the associated damage is treated as a crack layer (CL) 
(9). The damage parameter for simplicity is defined as an area of 
discontinuity surfaces per unit volume. A part of the CL within Which the 
rate of damage accumulation is positive, (^ > o) is defined as an active 
zone. The CL theory thus models fracture propagation as the active zone 
movements: translation and rotation as a rigid, isotropic expansion and 

distortion for a small active zone. There are potential energy release rates 
corresponding to each of these elementary movements. Energy release rates 
corresponding only to the active zone translation, the -if-integral, and 
expansion, the M-integral are considered here. The rate of the damage density 
resulting from these two movements of active zone is ilO,ll]. 

j> X* dkf> 

where v^ and e stands for the rates of translation in direction tangent to 
the crack path and expansion, respectively, and X is a position vector with 
respect to the center of isotropic expansion. 

Potential energy release rate P associated with damage growth is 
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where 7f is e potential energy density and integral is taken over the active 
zone. Upon substitution of (1) this relation takes the fora 


f> V*b#dv 


( 3 ) 


where k * 1,2 

The first integral in the right-hand side of (3) is a translational energy 
release rate (J-integral) , the second one is an expanslonal energy release 
rate (M- integral ) . 

2. Translational Energy Release Rate in Thermoelasticity 

The stress tensor £*and the absolute temperature T constitutes a 
conventional set of parameters for a thermos last ic medium. This set, extended 
by adding the damage parameter j} , is chosen as a system of thermodynamic 
state parameters 

W If] 


The corresponding thermodynamic potential is a Gibbs' free energy density 
which is usually designated in continuum mechanics as the potential energy 
density 


ir= 


(4) 


Then, the following constitutive equations for thermoelasticity are held: 
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(5) 
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h * CoUST 

(6) 


3 


where S is the entropy, end ^ is the totel strain tensor. According to Bq. 
(3). 


a-Sf^Jv 


Since (4) 






3T 


(7) 


( 8 ) 


Using local constitutive equations (5) and (6) this expression can be 
rewritten as 


d*P c d* f-fyfcj Vi,*) + S d*T C i,4,K - 4 a) 


(9) 


Hero, f is the strain energy density, or the specific Helmholtz free energy, 
and the component of the displacement vector. 

For flC being a linear coefficient of thermal expansion the specific 
entropy for thensoelasticity 


- oC Ot 


kk 


( 10 ) 


Substituting (9) and (10) into expression (7) for and using Gauss theorem 
we finally arrive at the integral 


F 1 ^ l f n l “ ^ij n j u i,l ld ^^ c ^ 0kk 0,1 


dv 


( 11 ) 
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Here, is j-th component of the unit normal to an integration path 
P^ 9 * * - T q , where T q is the absolute temperature for the initial 

equilibrium state of the system when the strain is everywhere zero. The first 
term in (11), i.e., integral along the path P repeats that for elastic case 
with the fields f , £* and V corresponding to the thermoelastic equilibrium. 

The second term, i.e., the area integral, cannot be in general converted into 
a line integral. However, such a conversion can be easily done for linearly 
varying temperature field. 

3. Expans ional Energy Release Rate 

Expression for the H-integral follow from the second part of the 
right-hand side of (3): 


M X* <Jv (1 

V 

Using the same approach like in the previous section the expression for this 
integral is obtained in the following form. 


M s ^\ [ X* X* P * 0J| ( 0 t 0,*) J V 

» 


(13) 


where i, j, k and Jl *1,2. 

The position vector is taken with respect to the center of isotropic 
expansion. 

The M-integral consists of two parts similarly to the J-integral (11): one 
is a line integral and another is an area integral. 

4. Discussi on 

Since p is not zero within the active zone only the following integral 
over an arbitrary area V is reduced to 
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( 14 ) 


A#'*- S*'* 

It leads to the following relations 

v '* v* 

and 

S If Xk v ^ * S af ** ^ 

v' v5» 


These imply the invariancy of J- and H- integral with respect to the domain 
of integration V V^. As it was shown above the area integral can be 
decomposed into the linear and area integrals. The linear integral formally 
coincides with that for the elastic case, however, for thermoelasticity it is 
not path independent. Only the sum of the linear and area integrals in (11) 
and (13) obeys the conservation law. 

The expressions (11) and (13) agree with corresponding integrals of [12] 
and differ from all the rest mentioned above expressions for energy release 
rate. 
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